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AUTONOMIC  PHYSIOLOGICAL  DATA  ASSOCIATED 
WITH  SIMULATOR  DISCOMFORT 


INTRODUCTION 


There  has  existed  a  controversy  concerning  the  relationships  of 
autonomic  functions  to  no'^io^sickness  for  at  least  two  decades. 
Money  exhaustively  reviewed  the  signs  and  symptoms  of  motion 
sickness  (1) .  For  those  organs  innervated  by  the  autonomic 
nervous  system  <ANS) ,  he  noted  that  conflicting  observations 
(e.g.,  both  increases  and  decreases  in  heart  rate)  had  been 
reported . 

Money  argued  (1),  in  part,  that  since  (a)  visual*^e^ibular 
interactions  reside  in  the  central  nervous  system  (CNS) ,  (b) 
emesis  has  strong  somatic  muscular  components  (diaphragm  and 
abdominal  wall) ,  (c)  skin  pallor  can  be  induced  at  the 
sympathetic  ganglion  by  'acetylcholinesterase,  and  (d)  that  nausea 
may  reflect  reactions  to  CNS  events  alone,  then  one  may  not 
assume  that  motion  sickness  is  a  dysfunction  within  the  autonomic 
nervous  system  (ANS) .  An  alternative  h^^othesls  is  that  motion 
sickness  is  a  CNS  dysfunction  which  usually  has  sequelae 
involving  the  ANS. 


Subsequently,  Graybiel  and  Lackner  were  unable  to  find  systematic 
autonomic  indicators  of  motion  sickness  within  or  across  subjects 
(2) .  The  measures  available  to  the  investigators  were  somewhat 
non-specific:  heart  rate,  blood  pressure,  and  body  temperature 
are  all  subject  to  many  intrinsic  and  extrinsic  influences. 
However,  Cowings  et  al.  found  systematic  autonomic  responses, 
then  suggested  that  stability  of  ANS  responding  exists  within 
subjects  (3,  4).  They  used  somewhat  more  specific  measures  than 
Graybiel  and  Lackner,  including  skin  conductance  and  finger  pulse 
volume,  among  others. 
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In  the  neantlne,  Stem  and  his  colleagues  at  Pennsylvania  State 
University  showed  across-subject  consistency  in  the  hypergastrlc 
response  to  motion,  and  anecdotal  evidence  of  a  correlation 
between  hypergastrla  and  subjective  reports  of  gastric  discomfort 
and  nausea  (5,  6,  7) . 


This  somewhat  confused  state  of  affairs  concerning  ANS  responses 
to  motion  sickness  presented  a  problem  to  applied  researchers 
associated  with  the  U.S.  Army's  fixed-base,  virtual-world 
helicopter  simulator  at  NASA  Ames  Research  Center.  In  the  realm 
of  engineering,  the  device  provided  a  providential  proving  ground 
for  controls  and  displays  research.  Unfortunately,  a  simulator 
slclcness  rate  of  20  to  40%  was  expected  for  experienced  pilots 
operating  the  device  (8) . 

We  report  here  the  development  of  a  physiological  monitoring 
cap2U3lllty  for  the  Army's  Crew  station  Research  and  Development 
Facility  (CSROF)  and  some  preliminary  physiological  data.  The 
objective  this  effort  was  to  demonstrate  sensitivity  of  the 
physiological  measures  In  this  simulator  to  self-reported 
simulator  sickness. 


METHOD 


Eight  channels  of  data  were  gathered  through  a  i2-blt  (4096-step} 
analog-to-digital  system  (Qua tec,  Akron  OH,  Model  SAC-12)  with  a 
total  sampling  rate  of  33  mHz.  Six  channels  were  autonomic 
physiological  data:  electrocardiogram  (ECG) ,  forearm  skin 
conductance  level  (SCL) ,  electrogastrogram  (EGG) ,  chest 
clrctUBference  for  ventilatory  rate  (ty) ,  fingertip  pulse  (BVP), 
and  fingertip  skin  temperature  (T,,^) .  One  channel  was  an  event 
marker  for  flight  segment  time  registration.  The  last  channel 
was  unused. 
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signal  Acquisition 


The  six  physiological  signal  conditioners,  built  by  NASA-Anes, 
were  Isolated  at  the  power  supply  (voltage  converter)  level. 

They  used  front-end  band-pass  filtering  specific  to  the  signal  of 
interest,  and  amplified  each  signal  to  the  ±5  v  p-p  range.  The 
differential  amplifiers  used  for  ECG  and  EGG  had  common  mode 
rejection  ratios  exceeding  80  db. 

The  adhesive  spot  electrodes  used  for  ECG,  EGG,  and  SCL  were  of  a 
standard  Ag/AgCl  clinical  ECG  variety  with  snap  connectors 
(Vermont  Medical,  Inc.,  Bellows  Falls,  VT  05101,  Re-order  no. 
A-10005) .  These  electrodes  were  recessed,  with  a  thin  piece  of 
foam  between  the  skin  and  the  electrode  surface,  to  minimize 
motion  artifact.  Clinical  electrolyte  (TECA  Corp. , 

Pleasantville,  NY  10570,  Cat.  no.  822-201210,  NDC  0324-0045-06) 
was  used  at  the  electrode-skin  interface. 

Creation  of,  and  access  to,  the  digitized  data  files  was  mediated 
by  the  software,  CODAS  (Data-Q  Instruments,  Inc.,  Akron  OH, 
version  3.0),  which  ran  at  a  total  sampling  rate  of  800  Hz  (100 
Hz  per  channel) .  CODAS  also  supported  visual  data  display  and 
ASCII  or  Lotus  123  file  creation. 

Electrogastrogram  -  The  EGG  signal  conditioner  analog  band  pass, 
0.02  to  0.25  Hz  (1.2  to  15  cpm) ,  allowed  the  detection  of 
"hypergastrla"  in  the  4  to  9  cpm  range  (6) .  We  used  the  two 
active  EGG  electrode  placements  of  stern  et  al.  (6):  (l)  lO  cm 

above  the  navel  and  6  cm  left  of  the  centerline,  and  (2)  4  cm 
edjove  the  navel  on  the  centerline,  with  a  reference  electrode 
placed  on  the  left  mid-axillary  line,  halfway  between  the  iliac 
crest  of  the  hip  and  the  bottom  of  the  rib  cage. 

Skin  Conductance  Level  -  Though  the  hypothalamic  control  of 
eccrine  thermoregulatory  sweating  is  well  described,  and  the 
responsiveness  of  the  SCL,  mediated  by  eccrine  sweat  gland 
activity,  to  motion  discomfort  has  been  noted  (3),  the  complete 
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picture  of  the  neural  circuitry  underlying  eccrine  gland  activity 
is  not  clear  (9) .  However,  decades  of  experience  suggest  that 
the  SCL  may  provide  reasonable  indications  of  emotion-triggered 
activity  in  the  sympathetic  branch  of  the  ANS. 

Due  to  the  rec[uirement  that  the  pilot  use  his  hands  to  fly,  and 
the  impractical ity  of  electrode  placements  on  the  sole  of  the 
foot,  the  two  SCL  electrodes  were  attached  to  the  skin  of  the 
medial  wrist.  They  were  placed  laterally  to  one  another,  just 
medial  to  the  anterior  annular  ligament.  The  skin  was  not 
abraded. 

The  SCL  signal  conditioner  placed  a  15  Hz,  constant  1.25  v  p-p 
square  wave  across  the  skin  and  a  100  kohm  current  limiting 
resistor  in  series  with  the  skin.  SCL  variations  were  viewed  in 
the  frequency  range  dc  to  2  hz  as  the  voltage  required  to  hold 
current  constant  across  the  combined  resistive  impedances.  The 
current  flow  through  the  series  resistances  varied  as  allowed  by 
skin  conductance,  the  Inverse  of  resistance.  The  signal 
conditioner  voltage  output  was  nonlinear,  and  was  corrected  by 
the  formula. 
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where  gain  =  1,453,682  and 

offset  *  -5.73354  volts 

The  SCL,  in  mhos  (a  unit  of  electrical  conductance) ,  was 
calculated  as  the  reciprocal  of  Conductance  has  been 

preferred  to  resistance  in  electrodermal  response  analysis,  due 
to  a  more  Gaussian  distribution,  for  a  number  of  years  (lO). 

Electrocardiogram  and  Vagal  Tone  -  The  ECG  provides  a  window 
into  electrical  conduction  characteristics  of  the  cardiac 
neuromuscular  structure.  It  is  also  one  technique  used  for  the 
measurement  of  heart  rate.  Both  conduction  and  rate  are  modified 
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by  autonomic  humoral  mechanisms,  and  rate  is  modified  by 
intrinsic  myocardial  length-tension  relationships  (the 
Frank-Star ling  mechanism) . 

In  addition,  rate  is  influenced  by  autonomic  parasympathetic 
fibers  travelling  in  the  tenth  cranial  nerve,  the  vagus. 

Efferent  cardiac  fibers  originate  in  the  medulla  and  travel  in 
the  vagus.  Also,  the  efferent  arms  of  reflex  circuits 
originating  from  pulmonary  and  cardiac  receptors  travel  in  the 
vagus.  Rate  is  also  Influenced  by  autonomic  sympathetic  fibers 
originating  in  the  cardiac  nuclei  of  the  sympathetic  trunk. 

The  reflex  effects  of  ventilatory  rate  and  depth  on  heart  rate 
are  well  documented  (11) .  Scientific  literature  published  over 
70  years  ago  showed  "(1)  that  respiratory  sinus  arrhythmia  is 
mediated  by  the  vagus;  (2)  that  the  amplitude  of  respiratory 
sinus  arrhythmia  is  related  to  the  functional  status  of  the 
cardiac  vagi  (l.e.,  cardiac  vagal  tone);  and  (3)  that  an 
individual  with  pronounced  respiratory  sinus  arrhythmia  has 
specific  behavioral  characteristics."  (12) 

Investigators  have  long  sought  a  window  on  vagal  activity  by 
quantifying  respiratory  sinus  arrhythmia.  Recently,  several 
laboratories  independently  identified  three  bandwidths  of  heart 
rate  variation  (13).  In  humans,  studies  have  indicated 
relationships  between  heart  rate  variations  around  0.10  Hz  and 
mental  work  load  (14,  15).  The  responses  of  the  .12  to  4.0  Hz 
component  of  sinus  arrhythmia  to  atropine  (16)  and  in  the 
presence  of  alpha-  and  beta-adrenoreceptor  blockers  indicates  its 
mediation  by  the  parasympathetic  branch  of  the  autonomic  nervous 
system,  in  particular,  the  vagus  nerve  within  which  the 
applicable  fibers  tun  (17). 

The  Vagal  Tone  Monitor  (VTM;  Oelta-Biometrlcs,  Bethesda,  MD) 
provided  automated,  on-line  data  reduction  from  the  raw  ECG  to 
one  CO  ponent  of  respiratory  sinus  arrhythmia.  We  used  the  VTM 
to  assess  vagal  activity,  or  vagal  tone  (VT) ,  monitoring  the 
middle  component  of  sinus  arrhythmia,  the  0.12  to  0.40  Hz  band. 
This  provided  us  with  an  index  of  parasympathetic  activity.  The 
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VTM  also  provided  a  cardiotachometer  function  from  which  we 
extracted  the  directly-measured  heart  period  (HP) ,  the  reciprocal 
of  heart  rate.  HP  was  selected  as  a  measure,  rather  than  heart 
rate,  for  its  more  Gaussian  distribution. 

The  ECG  signal  was  derived  from  electrodes  placed  to  generate  the 
"CR-S**  lead  (18) .  The  CR-5  lead  is  an  electrode  placement  which 
minimizes  movement  and  EMG  artifacts  in  the  ECG  and  produces  a 
lead-II-appearing  (large  R  wave)  tracing.  The  lead  was  described 
by  Blackburn,  who  investigated  the  characteristics  of  22 
different  bipolar  chest  leads  (19).  The  "CS"  lead  was  recorded 
from  the  manubrium  of  the  sternum  or  from  the  right  clavicle  to 
the  C5  position  (left  anterior  axillary  line  at  the  5th 
intercostal  space.  It  was  the  most  sensitive  lead  for  the 
detection  of  S-T  segment  changes  after  exercise.  We  used  the 
standard  electrocardiography  color-code  and  nomenclature  for 
electrode  connectors.  The  white  lead  (RA) ,  to  minimize  motion 
and  EMG  artifact  from  the  clavicular  portion  of  the  pectoralis 
major,  was  placed  when  possible  upon  the  manubrium  of  the 
sternum,  just  below  the  notch  at  the  top  of  the  sternum.  Since 
pectoral  activity  in  this  experiment  vas  mlnir.al,  positions  at 
the  level  of  the  manubrium,  but  displaced  to  the  subject's  right, 
were  used  when  chest  hair  patterns  were  a  problem. 

The  black  lead  (lA)  was  attached  directly  below  the  left  nipple 
in  the  5th  intercostal  space.  We  used  finger  pressure  on  the  rib 
cage  to  find  the  intercostal  space.  Placement  directly  over  a 
rib  might  have  caused  motion  artifact.  The  RA  and  LA  connections 
led  to  the  active  inputs  of  a  differential  amplifier  with  analog 
front-end  band-pass  filtering  at  0.5  to  100  Hz.  The  green  (RA; 
reference  electrode)  lead  was  attached  slightly  forward  of  the 
right  mid-axillary  line,  halfway  between  the  iliac  crest  of  the 
hip  and  the  bottom  of  the  rib  cage.  Placement  on  the  iliac  crest 
or  on  the  lower  ribs  might  have  caused  motion  artifact. 
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Subjective  and  Other  Evaluations 


The  pilots  were  asked  to  rate  their  motion  discomfort  (MD)  on  a 
7-polnt  scale,  on  which  1  was  "Normal,  symptom  free,"  7  was 
"Severe  discomfort,  I  am  unable  to  continue,"  and  the  values  2 
through  6  were  not  anchored  (8) .  The  ratings  were  collected  once 
every  five  minutes  while  the  pilots  were  in  the  simulator.  Also, 
the  pilots  completed  a  motion  history  questionnaire  and  a 
simulator  side  effects  questionnaire  (SSEQ) .  In  addition,  they 
received  two  postural  equilibrium  assessments:  a  stand-on-one- 
leg,  eyes-closed  test  (SOLEC)  and  a  walk-on-straight-line,  eyes- 
closed  test  (WOFEC) .  The  measures  of  interest  in  the  latter  were 
the  total  time  of  balance  across  three  trials  (SOLEC)  and  the 
total  number  of  heel-to-toe  steps  accomplished  prior  to  loss  of 
balance  in  one  trial  (WOFEC) . 


Data  Reduction 


All  channels  of  analog  physiological  data  were  sampled  at  100 
82unples/sec,  then  reduced  in  the  following  manner.  The  skin 
conductance  level  (SCL)  data  were  summarized  as  mean  values  for 
30-second  epochs.  Then  the  linear  trend  was  calculated,  using 
the  least  squares  method,  for  each  session,  and  subtracted  from 
the  raw  data.  The  residual  30-second  epoch  means  were  used  for 
further  analyses.  These  30-second  means  fell  into  the  domain  of 
tonic,  rather  than  phasic,  SCL  measures. 

Cardiac  interbeat  intervals,  or  heart  periods  (HP),  were  also 
summarized  as  mean  values  for  30-second  epochs.  The  variance  in 
cardiac  Interbeat  interval  was  partitioned  such  that  variance  in 
the  frequency  band,  0.12  to  0.40  Hz,  was  reported  each  30  seconds 
as  vagal  tone  (VT) .  The  VT  data  were  smoothed  by  averaging 
together  four  sequential  epochs  (two  minutes) ,  creating  a  moving 


7 


average  lagged  30  seconds  at  a  tine  and  labelled  by  the  fourth 
epoch  In  the  average. 

For  subsequent  analyses  involving  the  EGG  data,  the  HP,  SCL,  and 
VT  data  were  reduced  to  1-minute  epochs  by  averaging  together 
sequential  pairs  of  30-second  epochs,  without  overlap. 

The  digitized  EGG  data  for  one  minute  (6000  samples)  were  shifted 
to  zero  mean,  tapered  at  both  ends  (10%  taper)  (20) ,  and 
subjected  to  a  discrete  Fourier  transform  (MATLAB,  The  Math 
Works,  Inc.,  South  Natick,  MA) .  The  output  of  the  transform 
included  raw  energy  estimates  and  phase  estimates  in  each 
frequency  bin  from  1  through  13  cycles/minute.  The  energy 
estimates  were  reduced  to  an  index  of  normal  gastroenteric 
activity  (EGG3;  3  cycles/min)  and  of  hypergastria  (HG;  mean  of  4 
through  9  cycles/min)  (5) .  The  raw  power  data  were  transformed 
to  logarithms  to  provide  more  Gaussian  data  distributions.  The 
two  indices  were  smoothed  by  averaging  together  four  sequential 
minutes  of  log  power  data,  creating  a  moving  average  lagged  one 
minute  at  a  time  and  labelled  by  the  fourth  minute  in  the 
average . 

The  MD  ratings  were  generalized  across  5-minute  periods.  Each 
rating  value  was  assigned  to  the  minute  in  which  it  was  reported, 
and  to  the  two  minutes  immediately  preceding  and  immediately 
following  the  report.  For  example,  the  report  ”3”  in  minute  20 
and  "4”  in  minute  25  caused  minutes  18,  19,  22.,  21,  and  22  to  be 
assigned  the  value  "3,"  and  minutes  23,  24,  25.,  26,  and  27  to  be 
assigned  the  MD  rating  ”4."  When  a  2-rank  change  occurred,  the 
rating  value  was  assigned  to  the  minute  in  which  it  was  reported, 
and  the  one  minute  immediately  preceding  or  immediately  following 
the  report,  depending  on  the  sequence  of  change.  The  other,  more 
distant  minute  preceding  or  following  the  report  was  assigned  the 
next  report  level.  For  example,  the  report  ”2”  in  minute  20  and 
"4"  in  minute  25  caused  the  sequence  2-2-2,-2-2-3-4-4.-4-4 ,  for 
minutes  18  through  27,  respectively.  No  3-rank  changes  occurred 
across  sequential  5-min  periods. 
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Flight  Conditions 


The  pilots  flew  in  four  experimental  conditions,  created  from  two 
experimental  factors  with  two  levels  each.  The  conditions  were 
low-effort  maneuvering  at  100  ft  above  ground  level  (AGL) ,  low- 
effort  maneuvering  at  400  ft  AGL,  high-effort  maneuvering  at  100 
ft  AGL,  and  high-effort  maneuvering  at  400  ft  AGL.  The  flight 
conditions,  pilot  performance  data,  equilibrium  data,  and  motion 
discomfort  data  are  reported  fully  elsewhere  (21) .  The  full 
analysis  indicated  that  only  the  altitude  factor  had  a 
significant  effect  upon  motion  discomfort.  Thus,  a  subset  of 
data  were  selected  for  an  analysis  here  of  physiological  data 
which  would  reflect  the  effects  of  the  altitude  manipulation. 

The  data  came  from  two  independent  groups  of  subjects  with  n=6  at 
100  ft  AGL  and  n=7  at  400  ft  AGL. 


Statistical  Analyses 


Graphic  displays  of  the  individual  and  group  physiological  data 
were  used  to  search  for  apparent  differences  related  to  altitudes 
and  related  to  motion  discomfort  ratings.  For  the  group  data, 
when  a  MD  contrast  was  desired,  an  approximate  median  split 
approach  was  used.  The  subjects  were  separated  into  a  low 
ratings  group  (maximum  ratings  of  "I"  through  "3")  and  a  high 
ratings  group  (maximum  ratings  of  ''4''  through  ”7”)  . 

The  data  set  was  not  complete  enough  nor  extensive  enough  to 
attempt  to  partition  variance  among  MD,  subjects  and  error 
factors  using  an  analysis  of  variance.  However,  the  data 
contributin  j  to  the  graphic  views  of  the  data  allowed  a 
simplistic  sensitivity  assessment  of  the  different  physiological 
measures.  The  metric  used  for  sensitivity  assessment  was  the 
ratio  of  the  average,  between-condltlon  difference  (d)  and  the 
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grand  standard  deviation  (s) .  The  latter  was  calculated  from  all 
availedale  observations  of  the  variable. 

Discriminant  analyses  (StatGraphics,  STSC,  Inc.,  Rockville  MD, 
version  4.0)  were  used  to  search  for  predictive  relationships 
between  motion  discomfoirt  ratings  and  a  multivariate  set  of 
physiological  activities.  This  approach  was  used  for  individual 
and  for  group  data.  For  discriminant  analyses  using  group  data, 
5-ininute  MD  ratings  were  categorized  as  low  (ratings  of  "I”)  , 
medium  (ratings  of  ”2"  and  ,  and  high  (ratings  of  ''4”  through 
"T") .  For  discriminant  analyses  using  individual  data,  the  5- 
mlnute  MD  ratings  were  not  combined.  Five  physiological 
variables  were  used  to  predict  MD:  HP,  SCL,  HG,  EGG3  and  VT. 

Each  variable  was  represented  as  within-subject,  within-session, 
within-variable  standard  scores  (mean=0,  s=l) .  The  variable, 
time,  was  purposely  avoided  as  a  predictor  to  allow  a  full 
exploration  of  the  predictability  of  MD  from  these  physiological 
data.  Not  enough  data  were  available  to  both  "train”  and 
validate  the  discriminant  functions  within  or  across  subjects. 


RESULTS 


We  recovered  usable  SCL  data  from  ii  subjects,  HP  and  VT  data 
from  13  subjects,  EGG  data  from  11  subjects,  and  MD  ratings  from 
9  subjects,  as  shown  in  Table  1.  Eight  subjects  provided 
complete  sets  of  data. 

Standard  score  data  for  the  five  physiological  variables, 
compiled  across  the  eight  subjects  with  complete  data  sets,  were 
relatively  Gaussian  for  each  variable  (Figure  1) . 
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TABLE  1.  DATA  RECOVERY,  BY  SUBJECTS  AND  VARIABLES,  AND  MAXIMUM 
MOTION  DISCOMFORT  (MD)  RATINGS.  SCL,  SKIN  CONDUCTANCE  LEVEL;  HP, 
HEART  PERIOD;  VT,  VAGAL  TONE;  EGG,  ELECTROGASTROGRAM ;  ALL,  ALL 
DATA  AVAILABLE. 


Sub  j . 

Alt 

SCL 

HP  & 

EGG 

MO 

Max  MD 

All 

VT 

S02B 

HIGH 

X 

X 

X 

X 

1 

S03A 

HIGH 

X 

X 

X 

X 

4 

X 

S04A 

LOW 

X 

X 

X 

X 

2 

X 

SOSA 

LOW 

X 

X 

... 

.  •. 

S06A 

HIGH 

X 

X 

X 

..... 

... 

... 

S07A 

LOW 

X 

1 

1 

J 

... 

... 

... 

S08A 

HIGH 

X 

X 

X 

X 

6 

X 

SOSA 

LOW 

X 

X 

X 

X 

5 

X 

SlOA 

HIGH 

X 

X 

X 

2 

... 

S12A 

HIGH 

X 

X 

X 

X 

3 

X 

S14A 

HIGH 

X 

X 

X 

X 

3 

X 

SISA 

LOW 

X 

X 

X 

... 

... 

... 

S16A 

LOW 

X 

X 

X 

X 

4 

X 

Graphic  Views 


The  group  mean  HP  data  revealed  a  systematic  difference  for  the 
altitude  treatment  (Figure  2) .  The  heart  rate  of  the  low- 
altitude  group  (n-6)  was  consistently  higher  (shorter  HP)  than 
the  heart  rate  of  the  high-altitude  group  (n»7) .  The  mean  HP 
difference  between  the  low  and  high  altitude  conditions,  for 
epochs  1  through  60,  was  84.4  msec.  The  grand  standard  deviation 
of  all  obsezrvatiens  (1522  observations)  in  the  supporting  data 
set  was  88.5  msec.  The  ratio,  d/s,  was  .95. 

The  group  mean  HP  data  also  revealed  a  systematic  difference  for 
the  MD  median  split  (Figure  3).  The  heart  rate  of  the  low-rating 
group  (n«5)  was  consistently  lower  (longer  HP)  than  the  heart 
rate  of  the  high-rating  group  {n»4) .  The  mean  d  for  HP  between 
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Figure  l.<-  Frequency  distributions  of  across>subject 
observations,  using  within-subject  standard  scores,  of  normal 
gastric  activity  (EGG3;  log  power),  vagal  tone  (VT) ,  hypergastria 
(HG;  log  power) ,  tonic  skin  conductance  level  (SCL;  mhos) ,  and 
heart  period  (HP;  milliseconds) .  Numbers  of  subjects  and 
observations  available  in  Table  2. 
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Figur«  2.-  Mean  heart  period  data  from  low  (100  ft  AGL)  and  high 
(400  ft  AGL)  altitude  groups,  and  one  grand  standard  deviation 
(Table  2)  above  the  low  altitude  mean. 
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1 


Figure  3.-  Mean  heart  period  data  from  low  (1-3)  and  high  (4-7) 
notion  dlsconfort  (MD)  rating  groups,  and  one  grand  standard 
deviation  (Table  2)  bolow  the  low  MD  mean. 
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low  and  high  MD  categories,  for  epochs  1  through  100,  was  47.6 
nsec,  and  the  ratio,  d/s,  was  .54. 

The  group  mean  SCL  data  revealed  no  systematic  differences  for 
the  altitude  treatment  or  for  the  MD  median  split.  The  group 
mean  VT  data  revealed  a  systematic  difference  for  the  altitude 
treatment  (Figure  4),  but  none  for  the  MD  median  split.  The  VT 
of  the  low-altitude  group  (n"6)  was  consistently  less  than  the  VT 
of  the  high-altitude  group  (n>7) .  The  mean  plot  of  the  low- 
rating  subjects  overlaid  that  of  the  high-rating  subjects.  The 
mean  d  for  VT  between  low  and  high  altitudes,  for  epochs  l 
through  60,  was  .86  xmits.  The  grand  VT  s  (1482  observations) 
was  .98  units,  and  the  ratio,  d/s,  was  .88. 

The  group  mean  EGG3  data  revealed  no  systematic  differences  for 
the  altitude  or  MO  comparisons.  The  group  mean  HG  data -revealed 
a  systematic  difference  for  the  altitude  treatment  (Figure  5) . 

The  HG  of  the  low-altitude  group  (n"4)  increased  steadily  through 
minute  29,  while  the  HG  of  the  high  altitude  group  (n-7)  held 
steady.  For  minutes  5  through  10  plus  minutes  25  through  30,  the 
HG  mean  log  power  d  was  .57,  the  grand  s  (524  observations)  was 
.84,  and  the  ratio,  d/s,  was  .68. 

The  group  mean  HG  data  also  revealed  a  systematic  difference  for 
the  MD  median  split  (Figure  6) .  The  HG  of  the  low-rating  group 
(n-5)  diminished,  while  the  HG  of  the  high-rating  group  (n-4) 
increased.  For  minutes  30  through  50,  the  HG  mean  log  power,  d, 
was  .84,  and  the  ratio,  d/s,  was  1.0. 

These  calculations  are  summarized  in  Table  2.  The  rank  order  of 
sensitivity  was: 

1.  HG  for  MO  (d/ 8-1.0) 

2.  HP  for  altitude  (.98) 

3.  VT  for  altitude  (.88) 

4.  HG  for  altitude  (.68) 

5.  HP  for  MD  (.54) 


Discriminant  Analyses 


Individual  discriminant  analyses,  performed  on  the  data  of  each 
of  the  six  subjects  with  maximum  MO  ratings  greater  than  "2," 
revealed  an  Idiosyncratic  predictability  of  MD  ratings.  The  six 
analyses  produced  solutions  ranging  from  two  to  five  functions. 
The  first  functions  explained  64  to  93%  (median-82. 5%) ,  and  the 
second  functions  explained  4.5  to  38%  (median-14. 5% ) ,  of  the 
variance  in  the  subjects'  data.  All  first  functions  were  highly 


Figure  4.-  Mean  vagal  tone  data  from  low  (lOO  ft  AGL)  and  high 
(400  ft  ACL)  altitude  groups,  and  one  grand  standard  deviation 
(Table  2)  above  the  low  altitude  mean. 
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Figure  5.-  Mean  hypergaetrla  data  from  low  (lOO  ft  AGL)  and  high 
(400  ft  AGL)  altitude  groups,  and  one  grand  standard  deviation 
(Table  2)  above  the  low  altitude  mean  (minutes  5-10)  and  below 
the  low  altitude  mean  (minutes  25-30) . 
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Figur*  6.-  Mttan  hypargastrla  data  from  low  (l>3)  and  high  (4-7) 
notion  discomfort  (MO)  rating  groups,  and  one  grand  standard 
deviation  (Table  2)  above  the  low  MD  mean  (minutes  30-50) . 
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TABLE  2.  SUHMARY  OF  MEAN  DIFFERENCE  (d)  AND  GRAND  STANDARD 
DEVIATION  (B)  DATA  FOR  FIVE  PHYSIOLOGICAL  MEASURES  IN  ALTITUDE 
(100  vs.  400  ft  AGL)  AND  MOTION  DISCOMFORT  (1-3  VS.  4-7) 
COMPARISONS. 


EGG 

HG 

HP 

SCL 

VT 

loq  power 

msec 

mhos 

units 

Grand  s 

1.12 

.84 

88.5 

2.40x10- 

6 

.98 

n 

11 

11 

13 

11 

13 

No.  of 

Observations 

505 

524 

1522 

1217 

1482 

Motion  Discomfort: 

d 

.84 

47.6 

d/s 

1.0 

.54 

— 

Altitude : 

—  j  '  .57 

1 

84.4 

.86 

JZS _ i 

■  ■■  1  1 

~  1  .68 

.95 

— 

.88 

significant  (p<.001)  by  Chi^.  Three  second  functions  were 
significant  (p<.05). 

No  single  physiological  variable  dominated  the  top  of  the  order 
completely  (Table  3) .  The  Friedman  two-way  analysis  of  variance 
(22),  applied  to  the  column-sums  of  Table  3,  produced  a 
moderately  signlfic2«nt  result  (Chi^^  ■  7.07,  df*4,  p<.20).  SCL 
entered  first  into  the  first  function  three  times,  and  second 
once.  HG  entered  firwt  into  the  first  function  once,  and  second 
twice. 

However,  these  individual  discriminant  analyses  provided  better 
than  expected  (50%)  correct  cla-^sifications  of  wl thin-subject  MD 
reports  (Table  4).  The  correct  classifications  ranged  from  80% 
to  100%  across  the  MD  categories. 

These  data  and  the  data  of  two  more  subjects,  combined  across 
subjects  as  within-subject,  within-session,  within-variable 
standard  scores  (mean-O,  s>l) ,  produced  correct  classifications 
of  near  or  less  than  50%  across  MD  groups.  This  analysis  was 
based  upon  the  eight  subjects  with  complete  data  sets  (Table  1) 
and  the  grouping  of  MD  ratings  into  three  categories,  low, 
medium,  and  high.  The  correct  classifications  for  these  groups 
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were  38%  (62  minutes),  29%  (39  minutes),  and  54%  (33  minutes), 
respectively . 


TABLE  3.  RANK  ORDERS  OF  FIRST-FUNCTION,  STANDARDIZED 
COEFFICIENTS  FROM  DISCRIMINANT  ANALYSES  FOR  SIX  SUBJECTS, 
PREDICTING  MOTION  DISCOMFORT  (MO)  FROM  NORMAL  GASTROENTERIC 
ACTIVITY  (EGG3),  HYPERGASTRIA  (HG) ,  HEART  PERIOD  (HP),  SKIN 
CONDUCTANCE  LEVEL  (SCL) ,  AND  VAGAL  TONE  (VT) . 


EBiB9i 

EGG3 

HG 

HP 

SCL 

VT 

S12 

5 

2 

3 

1 

4 

S14 

2 

3 

1 

4 

5 

S16 

S 

4 

3 

1 

2 

S3 

3 

2 

4 

1 

5 

S8 

3 

1 

5 

2 

4 

S9 

1 

2 

4 

Stuns 

19 

14 

20 

12 

25 

An  examination  of  with in-subject  intercorrelation  matrices 
revealed  no  systematic  pattern.  First,  we  expected  to  find  a 
relationship  between  VT  and  HP.  However,  the  Pearson  r  values 
were  small  and  variable,  ranging  across  subjects  from  -.029  to 
.420  (mean,  via  Fisher  z  •  .122,  df  from  37  to  59).  Second,  we 
expected  to  find  a  relationship  between  VT  and  HG.  However,  the 
Pearson  r  values  ranged  from  -.429  to  .442.  Generally,  HG  and 
EGG3  were  highly  negatively  correlated. 


CONCLUSIONS 


The  graphic  views  of  the  group  mean  data  suggested  some  autonomic 
patterns.  For  the  motion  discomfort  (MD)  median  split,  higher  MO 
ratings  seemed  to  be  associated  with  higher  heart  rates  and 
higher  HG  activities  (Figures  3  and  6) .  These  responses  to  MO 
align  with  observations  by  Cowings  et  al.  (3)  and  by  stern  et  al. 
(5). 
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TABLE  4.  CORRECT  CLASSIFICATIONS,  PROM  INDIVIDUAL  DISCRIMINANT 
ANALYSES,  BY  SUBJECT  (n»6)  AND  MOTION  DISCOMFORT  RATING  (MD) . 
NUMBER  OF  MINUTES  CORRECTLY  CLASSIFIED  AND  THE  RESPECTIVE  PERCENT 
CORRECT. 


For  tlio  altitude  comparison,  heart  period  (HP)  and  vagal  tone 
(VT)  seemed  to  take  on  pilot  mental  workload  connotations.  Vagal 
tone  was  lower  and  heart  rate  was  higher,  suggesting  higher 
workload  at  the  lower  altitude  (100  ft  AGL;  Figures  2  and  4) . 

This  apparent  workload  effect  was  consistent  with  the  views  of 
pilots  that  flights  at  lower,  nap-of-the-earth  altitudes  produce 
higher  workloads  than  flights  at  higher  altitudes. 

However,  the  data  supported  an  alternative  hypothesis  explaining 
the  HP  and  VT  differences  at  the  two  altitudes.  The  lower 
altitude  seemed  to  produce  a  slight  hypergascrlc  disturbance 
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(Figure  5) ,  consistent  with  the  hypothesis  that  a  lower  nap-of- 
^e-earth  altitude  would  produce  such  an  effect,  and  the  higher 
Incidence  of  notion  discomfort  at  the  lower  altitude  (21) .  We 
cannot,  from  these  data,  determine  whether  wor)cload  or  MO 
affected  HP  and  VT. 

The  discriminant  analyses  suggested  (l)  that  no  single 
physiological  variable  dominates  In  the  prediction  of  NO  across 
subjects,  (2)  that  physiological  variables  mav  predict  MO  when 
restricted  to  wlthin-subject  comparisons,  and  (3)  that 
physiological  varledales  may  not  predict  MO  when  combined  across 
subjects.  There  are  four  caveats.  First,  the  capability  to 
perform  stepwise  discriminant  analyses  would  have  allowed 
hypothesis  testing  when  each  physiological  variable  entered  Into 
a  function  (23),  allowing  a  more  definitive  approach  to  the 
question  of  Idiosyncrasies.  This  capability  should  be  added  to 
the  analysis  approach.  Second,  the  small  number  of  observations 
per  Individual  subject  (range  39  to  61)  probably  Introduced 
sampling  errors. 

Third,  we  were  unable  to  validate  our  discriminant  functions 
within  or  across  subjects,  or  with  data  from  another 
Investigation.  This  validation  worlc  Is  needed.  Finally,  the 
assumption  that  a  single  MD  report  may  be  generalized  across 
several  minutes  obscured  some  unknown  proportion  of  the  true 
classifications  of  minutes.  MD  reporting  once  per  minute  should 
be  attempted  In  the  future. 

The  objective  of  this  effort  was  to  demonstrate  the  sensitivity 
of  physiological  measures  In  this  simulator  to  self-reported 
simulator  sickness.  We  were  successful.  However,  we  would  hope 
to  develop  the  metrics  to  somewhat  greater  levels  of  sensitivity 
In  subsequent  work.  Figures  3,  5  and  6  and  Tables  2  and  3 
suggested  that  HP,  HG,  and  SCL  may  be  more  sens it. We  to  simulator 
slcluiess  than  the  other  two  measures  used  here  (VT  and  EGG3). 
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FOOTNOTES 

^Dr.  Miller  served  as  a  consultant  to  Monterey  Technologies.  His 
firm  is  Miller  Ergonomics,  Lakeside,  California. 
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